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The impacts of the M9.0 Tohoku Earthquake on deep-sea environment were investigated 36 and 98 days 
after the event. The light transmission anomaly in the deep-sea water after 36 days became atypically greater 
(—35%) and more extensive (thickness ~1500 m) near the trench axis owing to the turbulent diffusion of 
fresh seafloor sediment, coordinated with potential seafloor displacement. In addition to the chemical influx 
associated with sediment diffusion, an influx of 13 C-enriched methane from the deep sub-seafloor reservoirs 
was estimated. This isotopically unusual methane influx was possibly triggered by the earthquake and its 
aftershocks that subsequently induced changes in the sub-seafloor hydrogeologic structures. The whole 
prokaryotic biomass and the development of specific phylotypes in the deep-sea microbial communities 
could rise and fall at 36 and 98 days, respectively, after the event. We may capture the snap shots of 
post-earthquake disturbance in deep-sea chemistry and microbial community responses. 



Tectonic events such as earthquakes often trigger the typical environmental disturbances that include numer- 
ous topographic failures and landslides on the land and coast, and sudden outflows as well as cessations and 
changes in the chemical composition of groundwater and hot springs 1 " 3 . In past, disturbance in deep-sea 
environment after a M5 level earthquake swarm was detected as the turbulent diffusion of sediments and a 
manganese concentration anomaly in a seafloor observatory 4,5 . However, the earthquake swarm was still on very 
local scale. These M5 earthquakes are not entirely comparable with the millennial M9 earthquake such as the 
Tohoku Earthquake occurring on March 11, 201 1 6 " 9 . To determine the immediate and existing conditions in the 
deep ocean following the M9.0 earthquake and its aftershocks is a top-priority objective of multidisciplinary 
science and a basis for the socially necessary restoration of the marine ecosystem and fisheries in Japan. To this 
end, the initial chemical and microbiological observations 10 " 16 in the deep-sea water environments were made on 
April 15 and June 16, 2011 at the Tohoku Earthquake epicentral region, the Japan Trench. 

Results 

At all the stations (Fig. 1) at the time of 36 days after the M9 event, an intense light transmission anomaly (LTA) 
was found in the deep-sea bottom water (Fig. 2). The intensity and extent of LTA were greater and thicker, 
respectively, at the stations near the trench axis. The maximum LTA intensity and thickness were —35% and 
— 1500 m, respectively. These values were very unusual in deep-sea environments, even for hydro thermally active 
regions 17 . Seismological investigations indicated that the largest vertical seafloor displacement occurred around 
the stations R and Nl 7 ' 8 . Indeed, the actual seafloor displacement (horizontally 50 m and vertically 7 m) was 
observed in the trench slope between the large normal fault (Fig. 1) and the trench axis 9 . However, turbidity in the 
deep-sea bottommost water was still observed at all the stations investigated (except for Stn. JKEO) even 98 days 
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Figure 1 | [Main body] Map of study area. Sampling stations in this study (filled circles), past drilling sites (filled squares), a multi-channel seismic reflection 
survey line (bold line), and a hypocenter of the 3.11 M9.0 Tohoku Earthquake (filled star) are shown. [Insert] Potential sub-seafloor structure and tectonic 
properties of the landward slope of the Japan Trench around the epicentral region of the 3.11 Tohoku Earthquake (MY102). The sub-seafloor structure and 
tectonic properties are based on a multi-channel seismic reflection image 7 . Locations of CTD-hydrocast stations (vertical dash lines) are shown. 



after the M9 event (data not shown), although the turbidity anomaly 
pattern was generally similar in both the 36- and 98-days-after sam- 
ples. Thus, it seems very likely that the highly turbid and sizeable 
deep-sea water mass was produced by the turbulent diffusion of 
sediments during the seafloor ruptures and landslides, while it was 
hardly concluded whether the diffusion of sediments was caused by 
the M9 earthquake or by the frequently continuing aftershocks. 

As previously observed as a post- earthquake chemical disturbance 
in seawater associated with the diffusion of sediments 4 ' 5 , anomalous 
increase in dissolved manganese concentration relative to the back- 
ground level (<1 nmol/kg) was found in the 36-days-after deep-sea 
water samples showing significant LTA (Fig. 2). Although most of 
the dissolved manganese concentrations seemed to be relevant to the 
LTA values, several water samples showed highly variable dissolved 
manganese concentrations (Fig. 2). 



Methane is a frequently-used, effective chemical tracer for the sub- 
seafloor fluid inputs 18-20 whereas the post-earthquake methane influx 
from sub-seafloor reservoirs into deep-sea water has rarely been 
discussed 18 . In fact, the methane concentrations in the 36-days-after 
deep-sea bottom water were significantly higher than the usual back- 
ground level in deep-sea water (—0.5 nmol/kg) (Fig. 2). Because 
of the entirely oxic conditions indicated by dissolved oxygen (DO) 
concentration data (Fig. SI), it is very unlikely that the methane 
anomaly was induced by in situ microbial production (methanogen- 
esis). 

The stable carbon isotope composition (5 13 C) of deep-sea water 
methane was nearly constant at each station but differed among the 
stations in the 36-days-later samples (Fig. 2). Although it was known 
that microbial methanotrophic function increased the 5 13 C values of 
deep-sea water methane 21 ' 22 , a Keeling plot analysis 23 demonstrated 




LTA (%) Mn [nmol/kg] CH, [nmol/kg] 6 13 C (V) [nmol/kg] cell density (x10 4 ) [/ml] 

Figure 2 | Vertical profiles of light transmission anomaly, manganese and methane concentrations, stable carbon isotopic composition of methane, H 2 
concentration and microbial cell density in the 36-days-after deep-sea bottom water of all the stations. The altitude of sampling point above the seafloor 
is used as a vertical axis. LTA in the 36-days-after deep-sea water of each five station (R, Nl, N2, N3, and JKEO) are distinguished by colors (see legend). 
An open star in the cell density profile represents the average microbial cell density (1.1±0.12X10 4 cells/ml) in the deep-sea bottom water samples of the 
Japan Trench approximately 70 days after the M9 event (see SI for details). Open diamonds and triangle respectively represent samples taken at Stns. R 
and N2 98 days after the M9 event. 
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that most of the 5 13 C values of the deep-sea water methane could be 
interpreted as consequences of the simple dilution between each of 
two different types of endmember methane and ambient seawater 
methane (Fig. S2). One of the estimated endmember methane 
sources was highly 13 C- depleted and was obtained from Stn. R 
(— 80.3%o). The others were relatively 13 C-enriched and were ob- 
tained from Stns. Nl (-58.6%o), N2 (-56.6%o) and N3 (-65.5%o). 
The 5 13 C values of deep-sea water methane at the shallowest station 
of N3 demonstrated a considerably different mixing pattern from the 
other stations (Fig. S2). This result can be due to the potential mixing 
at this station between the sub-seafloor source of methane and the 
dissolved methane in the relatively shallower seawater, in which the 
13 C-enriched methane has been sometime found 24 . 

At fault movements, molecular hydrogen (H 2 ) can be mechano- 
chemically generated 25 and spread from the sub-seafloor fault inter- 
faces 26 . Prior to the investigation, we had hypothesized that trace 
amounts of the earthquake-generated H 2 might be detected in 
deep-sea water soon after the gigantic earthquake and its aftershocks. 
However, the deep-sea water H 2 was not abundant for detection 
(Fig. 2). In addition, we did not find deep-sea water samples contain- 
ing detectable sulfide (data not shown). Since it has been known that 
the H 2 and sulfide inputs in the deep seawater rapidly disappear 
through chemical and microbial processes 27 ' 28 , our data collected 
36 days after the M9 earthquake do not necessarily decline the pos- 
sible earthquake-induced influxes of these substances from the sub- 
seafloor environments. 

The possible responses of deep-sea planktonic microbial com- 
munities to the earthquake-induced environmental disturbances 
were also investigated at two periods of time after the M9 event. 
The microbial cell densities at several depths of the 36-days- after 
deep-sea water increased toward the seafloor were more than 
2.0X10 4 cells/ml (Fig. 2 and Table SI). The cell densities in the 



deep-sea bottom water samples (water depths of 5720 m and 
5000 m at Stn. R and 2940 m at Stn. N2) collected 98 days after 
the M9 event were 8.3X 10 3 , 9.9X 10 3 and 1.9X 10 4 cells/ml, respect- 
ively (Fig. 2 and Table SI). The cell densities in the 36 -days -after 
deep-sea bottom water samples were higher than those in the 98- 
days-after deep-sea bottom water samples at the similar depths 
(Fig. 2). The spatial difference in the microbial cell densities was also 
found in the comparison between the 36-days-after deep-sea bottom 
water samples at the trench regions (Stns. R, Nl, N2 and N3) and at 
the abyssal plain (Stn. JKEO) (Fig. 2 and Table SI). The 16S rRNA 
gene numbers in the whole microbial DNA assemblages showed 
quite similar patterns in spatial and temporal variation of microbial 
populations in the deep-sea bottom water environments as observed 
in the microbial cell densities (Table SI). In addition, when the 
proportion of archaeal 16S rRNA gene number in the whole prokar- 
yotic 16S rRNA gene number was compared among the deep-sea 
bottom water samples, the bottommost water samples at Stn. R at the 
times of 36 and 98 days after the M9 event exclusively showed the 
significantly lower proportion of archaeal 16S rRNA gene number in 
the whole prokaryotic 16S rRNA gene number (33.7 and 40.2%, 
respectively) than any other deep-sea water samples (73.4% on an 
average) (Table SI). 

As other responses of microbial community in the deep-sea water, 
spatial and temporal variability in the 16S rRNA gene phylotype 
composition was also examined in the deep-sea bottom water 
(Fig. 3). The whole prokaryotic phylotype compositions were 
broadly similar among the stations including Stn. JKEO, the various 
depths and the different sampling periods (Fig. 3). However, several 
signature phylogenetic groups were identified only in the deep-sea 
bottom water samples of the stations in the Japan Trench landward 
slope that were placed in the potentially largest seafloor displacement 
region: SUP05 phylogroup 34,35 within Gammaproteobacteria (Stns. 
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Figure 3 | Whole prokaryotic 16S rRNA gene phylotype compositions in the deep-sea water samples at all stations. The number in parenthesis shown in 
the right of the bar chart represents the number of 16S rRNA gene clones sequenced. Black reverse triangles indicate signature phylogenetic groups in the 
post-earthquake deep-sea water (see SI). Most of the names and abbreviations of the phylogenetic groups indicated at the right follow the Hugenholtz 
database 43 . The classification of Marine Group I crenarchaeota subgroups and other archaeal phylotypes follow a reference 44 . 
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Nl, N2, N3 and R); zeta-proteobacterial phylogroup 36 ' 37 (Stns. Nl 
and R); Desulfobacterales phylogroup 38 within Deltaproteobacteria 
(Stns. Nl and R), Desulfuromonadales phylogroup 39 within Delta- 
proteobacteria (Stns. R, N2 and N3), Arcobacter phylogroup 40 ' 41 
within Epsilonproteobacteria (Stns. Nl and R), DHVE5/Rice cluster 
V phylogroup 42 (Stns. Nl and R) and DHVE6 phylogroup 42 of the 
Archaea (Stns. R, Nl and N2). In addition, as demonstrated by the 
quantitative PCR estimation, the whole prokaryotic phylotype com- 
positions in the bottommost deep-sea water samples (both 36- and 
98-days-after samples) at Stn. R revealed the lower abundance 
of archaeal 16S rRNA gene phylotype populations than the other 
samples (Fig. 3). 

To understand the differences in the phylogenetic context of the 
post- earthquake deep-sea microbial communities, a statistical ana- 
lysis (UniFrac) 16 was conducted with the 16S rRNA gene phylotype 
compositions in the deep-sea water samples at the time of 36 and 98 
days after the M9 earthquake (Fig. S3). In the 36-days-after deep-sea 
water samples, the microbial phylotype compositions were broadly 
classified into 6 types: the bottommost water types of Stn. JKEO, Stn. 
R, Stn. Nl and Stn. N2, respectively, about 5000 m deep water type of 
Stn. R, and the overlying deep-sea water type at all the stations 
including the bottommost water sample of Stn. N3 (Fig. S3). When 
the microbial phylotype compositions were characterized in several 
98-days-after deep-sea water samples at the similar depth locations, 
the principal coordinates analysis (PCoA) revealed the significant 
time-course changes in the microbial phylotype composition (Fig. 
S3). The most drastic change was found in the phylotype composi- 
tion at the bottommost water of Stn. R and the phylotype composi- 
tion became much closer to the composition in the bottommost 
water at Stn. JKEO, a reference composition (Fig. S3). 

Discussion 

Potential chemical signs of the post-earthquake environmental dis- 
turbances were identified in the LTA anomaly, the concentrations of 
manganese and methane and the stable carbon isotope composition 
of methane in the deep-sea water environments 36 days after the M9 
event. The relationships between LTA intensity and manganese or 
methane concentration in deep-sea water were examined to estimate 
the input sources of manganese and methane 29,30 . The turbulent 
diffusion of sediments should accompany diffusive input of sediment 
pore-water chemistry into the overlying deep-sea water. Therefore, 
the manganese and methane concentrations would be linearly cor- 
related with the LTA intensity caused by the turbulent diffusion of 
shallow sediments. A reliable linear correlation of the LTA intensity 
with the manganese concentration was found in the deep-sea water 
samples at Stn. N2 and Stn. R when the extremely high-LTA (>20%) 
samples were excluded (Fig. 4a). The highly variable manganese 
concentrations in the high-LTA deep-sea water samples at Stn. R 
were probably due to the occurrence of complex leaching and 
scavenging processes of manganese with numerous numbers of 
suspended particles under the oxic condition. A similar linear cor- 
relation was identified only between the LTA intensity and the meth- 
ane concentration in the deep-sea water samples at Stn. R (Fig. 4b). 
The correlations strongly suggested that the manganese and methane 
in the deep-sea water samples at these stations would originate from 
the turbulent diffusion of shallow sediments. In contrast, the meth- 
ane (Stns. Nl and N2) and manganese concentrations (Stn. Nl) in 
the deep-sea water samples of several stations were quite deviated 
from the linear correlations (Fig. 4). These results indicated a possible 
input of these components from the sources other than the shallow 
sediments. 

Depth profile of stable carbon isotope composition of reserved 
methane in the sub-seafloor sediments was characterized at the sites 
of the Japan Trench landward slope near the epicenter of M9 earth- 
quake during Expedition 186 of the Ocean Drilling Project (ODP) 31 
(Fig. 1). The profile indicated the highly 13 C-depleted methane 
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Figure 4 | Relationship between (a) manganese and (b) methane 
concentrations and LTA value. ASW means ambient seawater. The blue 
and red broken lines for each panel are least- squares fitted lines for data 
from Stns. N2 and R, respectively. 



( 80%o) in shallow sediments and the increases of the 5 13 C values 

of subseafloor methane with depth 31 (Fig. 5). Similar vertical trends 
of the 5 13 C(CH 4 ) values were known also in other subseafloor meth- 
ane field 32,33 . The relatively 13 C-enriched methane, with a value of 

60%o, were found only in the sediments at depths of >800 m 

below the seafloor. Thus, potentially two types of endmember meth- 
ane inferred from the 5 13 C values would be derived from the different 
sedimentary reservoirs at different depths (Fig. 5). One source is the 
methane sorbed in the shallow fresh sediments, and is predominant 
in the deep-sea water at Stn. R. The other is the deep sub-seafloor 
methane introduced into the deep-sea water at Stns. Nl, N2 and N3. 
Although it is not yet known how such deep sub-seafloor methane 
migrate to the deep-sea water, Stns. N1-N3 are located in a more 
landward area of the Japan Trench than Stn. R, where the sub- 
seafloor environment is characterized by the tectonic extension 
and many normal faults 7 (Fig. 1). The massive earthquakes possibly 
induced hydrological changes by rapidly changing the pore-fluid 
pressure and local rock permeability, thereby altering the sub- 
seafloor fluid flows 1 ' 2 . Thus, it seems likely that the deep subsea- 
floor- derived methane would seep into the deep-sea water through 
earthquake-induced fluid pumping (Fig. 5). 
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Figure 5 | Vertical profiles of carbon isotope composition of methane in 
deep-sea bottom water in this study (altitude in meter [m]) and in sub-seafloor 
reservoir obtained from the expedition 186 of the Ocean Drilling Project 31 
(depth in meter below seafloor [mbsf]). As discussed in the literature, we can 
recognize distinct possible methane sources of turbulent diffusion of surface 
fresh sediment at Stn. R (a red arrow) and earthquake-induced pumping of 
deep sub-seafloor components at Stns. Nl, N2, and N3 (a blue arrow). 

Both the microbial cell densities and the 16S rRNA gene numbers 
in the deep-sea bottom water at Stns. R, Nl, N2 and N3 were evi- 
dently increased at the time of 36 days after the M9 event but 
decreased 98 days later to the levels that were comparable to the data 
at the reference station of JKEO (Fig. 2 and Table SI). These results 
strongly suggest that the increased microbial populations in the bot- 
tom water of the Japan Trench landward slope represent a response 
of the microbial community to the earthquake-induced envir- 
onmental disturbances. Although both the 36- and 98-days-after 
deep-sea bottom water samples were highly turbid and were likely 
influenced by the diffusion of sediments, the response of microbial 
biomass was found only in the 36-days-after deep-sea water envir- 
onments. This may imply that the earthquake-induced chemical 
impact in the deep-sea water such as discharges of subseafloor 
reduced chemical substances (e.g., CH 4 , H 2 , sulfide, ferrous iron 
and organics) could remain significant by the time of 36 days after 
the gigantic earthquake. 

The post- earthquake, spatial and temporal variation of the plank- 
tonic microbial communities in the deep-sea bottom water was also 
delineated by the 16S rRNA gene phylotype analysis (Fig. 3 and Fig. 
S3). Although the whole prokaryotic phylotype compositions were 
broadly similar among the various depths and the different sampling 
periods of all the stations, several signature phylogenetic groups 
were identified only in the deep-sea bottom water samples of the 
stations in the Japan Trench landward slope (Fig. 3). In addition, 
most of these signature phylogenetic groups were detected in the 



36-days-after deep-sea bottom water but became the undetected 
components at the similar deep-sea water environments 98 days after 
the M9 event (Fig. 3). Based on the phylogenetic affiliation, the 
signature phylogroups of the Bacteria and Archaea were related with 
potential chemolithoautotrophic or chemolithotrophic populations 
that thrive in the planktonic and sedimentary oxic-anoxic interface 
and the reducing habitats of the deep sea (see Supplementary Infor- 
mation). Thus, these signature phylotypes are probably derived from 
two sources: (i) indigenous benthic microbial components that 
inhabited the shallow sediments before the earthquakes and spread 
into the deep-sea water through earthquake-induced sediment dif- 
fusion and (ii) fresh planktonic components that dominated the 
microbial communities by responding to the earthquake-induced 
influxes of reduced energy sources in the deep-sea water. 

The post-earthquake spatial variation of deep-sea planktonic 
microbial community was illustrated in the comparison between 
the 36-days-after deep-sea bottommost water samples of all the sta- 
tions (Fig. S3). Although Stns. R, Nl, N2 and N3 were geographically 
located close together on the landward slope of the Japan Trench, 
the microbial phylotype compositions were considerably different 
among these stations, and the composition at the bottommost water 
of Stn. R (5800 m deep) was the most distantly related with those in 
the bottommost water of Stn. JKEO and in any other deep-sea water 
samples of the stations (Fig. S3). The distinctive microbial phylotype 
composition at the 36-days-after bottommost water of Stn. R is 
probably derived not only from the high frequency of emerging 
signature phylogenetic groups and but also from the relative enrich- 
ment of bacterial phylotype populations in the phylotype composi- 
tion (Table SI, Fig. 3 and Fig. S3). The microbial community in the 
bottommost water at Stn. JKEO is considered to represent a reference 
microbial community for the abyssal plain of the Pacific plate and to 
show little influence from the earthquake-induced environmental 
disturbances. Thus, the difference in the deep-sea microbial com- 
munity between Stn. JKEO and the other stations may represent the 
earthquake-induced environmental impact in the deep-sea bottom 
water habitats between the on-trench and off-trench locations, and 
even among the on-trench locations. The post-earthquake, temporal 
variation of deep-sea planktonic microbial community was inferred 
from the comparison between the 36- and 98-days-after deep-sea 
bottom water samples of Stns. R and N2 (Fig. S3). The phylotype 
compositions at the 98-days-after bottom water of these stations 
shifted closer to the composition of the 36-days-after bottommost 
water of JKEO (Fig. S3). The temporal shifts suggest that the micro- 
bial phylotype compositions in the 36-days-after deep-sea bottom 
water environments represent certain short-lived responses specif- 
ically to the chemical disturbance of the habitats induced by the 
gigantic earthquake. In addition, as observed in the variation of the 
microbial biomass and the signature phylogenetic components, the 
deep-sea bottom planktonic microbial communities may retrace the 
ordinal states relatively soon after the mitigation of post-earthquake 
environmental disturbance. 

Our results illuminated the transient and present natures of deep- 
sea chemical environments and microbiological communities after 
the 3.11 Tohoku Earthquake. The impacts of the M9.0 Tohoku 
Earthquake on deep-sea environment we detected were less severe 
than expected from the enormous energy release. However, the evid- 
ence of impact was spatially and temporally extensive. This evidence 
served to identify the earthquake-induced chemical influxes in the 
deep-sea water resulting from sub-seafloor fluid migration. A few 
oceanographic reconnaissances cannot reveal the detailed processes 
involved in earthquake-induced chemical disturbances and the cor- 
responding transition in the microbial ecosystem. Ongoing oceano- 
graphic expeditions and further multidisciplinary investigations will 
clarify the full impact of the massive earthquakes initiated from the 
March 2011 Tohoku Earthquake on the marine environments and 
ecosystems of northeastern Japan in the future. 
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Methods 

The Japan Trench is located on the northwestern margin of the Pacific Plate that is 
subducting beneath the northeastern Japan Arc (Supplementary Fig. SI). Multi- 
channel seismic reflection data at the epicentral region were acquired by the R/V 
Kairei (Japan Agency for Marine-Earth Science and Technology: JAMSTEC) in 1999 
(Line MY102 of KR99-08, analyzed in the reference 7 ) (Fig. 1). Vertical hydrocasts of 
the CTD-CMS (Conductivity Temperature Depth profiler with Carousel Multiple 
Sampling system) were conducted 36 days after the M9.0 event during the MR1 1-03 
cruise of R/V Mirai (JAMSTEC) at the epicentral region (Fig. 1). Stations Nl 
(38°10.6'N, 143°33.0'E, depth: 3502 m),N2 (38°08.7'N, 143°19.0'E, depth: 2954 m), 
and N3 (38°06.8'N, 143°05.0'E depth: 1948 m) are located on the western direction 
of the landward slope of a ridge associated with displacement along an outstanding 
normal fault that was considered to be the potentially largest slip of the Tohoku 
Earthquake 7 . Station R (38°12.5'N, 143°47.2'E, depth: 5715 m) is above a branch 
reverse fault ( Fig. 1 ) . The stations N2 and R were revisited during the YK 1 1 - E04 cruise 
of R/V Yokosuka (JAMSTEC) with another CTD-CMS. The station JKEO (38°00'N, 
146°30'E, depth: 5381 m) is the site of JAMSTEC Kuroshio Extension Observatory 
located in the abyssal plain of the Pacific plate. The Ocean Drilling Project (ODP) 
Sites 1150 (39°11'N, 143°20'E) and 1151 (38°45'N, 143°20'E) were the drilled sites 
during the expedition conducted in 1998 and the pore- water chemistry was already 
reported 31 . 

The wired CTD-CMS system of the MRU -03 cruise consisted of a CTD (SBE9 
Plus, Sea-Bird Electronics), a CMS (SBE32, Sea-Bird Electronics), 36 Niskin-X bottles 
(12-liter type, General Oceanics), a dissolved oxygen sensor (RINKO-III, JFE 
Advantech), and a light transmissometer (C-star 25-cm light-path type, WET Lab). 
The Light Transmission Anomaly (LTA), calculated from the difference between the 
in- situ light transmission value (Tr: %) and the value of the transparent layer at 
intermediate depth for each hydrocast, is used to describe deep-sea water turbidity. 
The CTD-CMS system of the YK11-E04 cruise consisted of a CTD (SBE11 Plus, 
Sea-Bird Electronics), a CMS (SBE32, Sea-Bird Electronics), and 12 Niskin-X bottles 
(12-liter type, General Oceanics). 

The seawater samples taken by the Niskin bottles were immediately subsampled 
into several optimized bottles for various geochemical and microbiological 
analyses 10 " 12 . For the analysis of manganese concentration 10 , the subsampled seawater 
in an acid-washed plastic bottle was filtered using a 0.22- um pore-size PTFE filter and 
acidified with nitric acid (TAMA Chemical) and analyzed by the luminol-H 2 0 2 
chemiluminescence detection method. For the analysis of methane 11 , sample seawater 
was subsampled into a 120 ml glass vial capped by a butyl-rubber septum after the 
addition of 0.5 ml HgCl 2 - saturated solution for poisoning. The methane concen- 
tration and carbon isotope ratio were simultaneously determined with a combination 
of purge and trap techniques and continuous-flow isotope ratio mass spectrometry. 
The stable carbon isotope ratio is presented in general delta notation on a permillage 
scale with respect to the Vienna PDB. For the analysis of molecular hydrogen con- 
centration 12 , the subsampling was conducted in a manner similar to that used for 
methane, except that a Teflon-coated butyl-gum septum was used rather than the 
untreated butyl-gum septum used for methane. The molecular hydrogen concen- 
tration was analyzed onboard using the headspace method with a gas chromatograph 
equipped with a trace-reduced gas detector (TRD-1: Round Science Inc., Japan) 12 
within six hours after the subsampling to avoid sample alteration during storage. If a 
hydrogen sulfide-like smell was detected in a sample, the seawater was subsampled 
and placed in a Nalgene bottle to quantify the H 2 S concentration using methylene 
blue colorimetry. 

Microbial cell counts were determined by a DAPI-staining direct count 
(Supporting Information). For extraction of microbial DNA, a portion (3 L) of deep- 
sea water was filtered with a 0.22-um-pore-size cellulose acetate filter (Advantec, 
Tokyo, Japan) and was preserved onboard at — 80°C. DNA was extracted by using the 
Ultra Clean Mega Soil DNA Isolation kit (MO Bio Laboratory, Solana Beach, CA, 
USA). Quantitative PCR of archaeal and entire prokaryotic 16S rRNA genes was 
performed using 7500 Real Time PCR System 1314 . Prokaryotic 16S rRNA gene clone 
analysis was conducted with another PCR experiment 15 (also see SI for details). To 
assess difference in phylogenetic context of the post-earthquake deep-sea microbial 
communities in the bottommost deep-sea water, an online tool, UniFrac 16 , was used 
for the principal coordinates analysis (PCoA). 
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